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Abstract 

Exposure to ozone (0 3 ) causes reduction both in the growth and yield of rice (Oriza sativa L). Commonly used 
Chlorophyll fluorescent measurements are not sensitive enough for short term exposure of 0 3 aiming an immediate 
assessments. Such a conventional method typically needs exposure over a few days to detect the influence. As an 
alternative method, we proposed a novel non-invasive, robust, real-time, optical Statistical Interferometric Technique 
(SIT to measure growth at an accuracy of 0.1 nm with a commonly consumed Japanese rice cultivar, Koshihikari. In 
the present study, we have conducted a repetitive 0 3 exposure experiment for three days under three different 
concentrations of 0 nl f 1 (control), 120 nl \'\ and 240 nl \'\ to investigate the damage and recovery strengths. As a 
measure to assess the effect and recovery from three consecutive day exposures of 0 3 , we measured the elongation 
rate (nm mm" 1 sec" 1 ) every 5.5 sec for 7 hours, and it revealed nanometric elongation rate fluctuations or Nanometric 
Intrinsic Fluctuations (NIF). Comparing the standard deviation (SD) of normalized nanometric intrinsic fluctuations 
(NNIF), which was normalized by that before the exposure, we found that drastic reductions under both 120 nl I" 1 and 
240 nl I" 1 0 3 concentrations. Reduction percentages were large under high 0 3 concentration of 240 nl I" 1 indicating the 
possibility of irreversible effect. However exposure to 120 nl I" 1 of 0 3 showed recovery on the 2 nd and 3 rd days. While 
SIT did reveal immediate effect based on an observation for a few hours, the visible foliar effect could be observed only 
after a week. Hence, the technique could provide a way for fast assessment of effect and recovery due to cumulative 
exposure of 0 3 and hence the tolerance as well as the vitality of plant. 

Keywords: Crop; Rice (Oriza sativa L); Growth; Plant growth measurement; Atmospheric pollution; Optical 
interferometric technique; Ozone; Chlorophyll fluorescence; Environmental effects; Instantaneous growth; Yield 



Introduction 

Tropospheric ozone (0 3 ) is an air pollutant generated by 
the photo oxidation of precursor gases, such as NOx 
and volatile organic compound (VOC) (Ainsworth 2008; 
Ashmore 2005; Darrall 1989). Further, 0 3 concentrations 
exceeding 80 nl T occur regularly, even in rural areas of 
India, China, and Japan (Beig et al. 2007; Wang et al. 
2007; Kobayashi et al. 1995; Kobayashi and Okada 1995). 
According to projections, 0 3 pollution may cause rice 
yield losses of up to 16% with no change in agricultural 
practices (Ainsworth 2008), which would put food security 
in Asia at a substantial risk. 
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It is known that 0 3 causes biochemical and physio- 
logical changes leading to the inhibition of photosynthesis 
and a consequent decrease in plant growth. 0 3 enters the 
leaf through stomata, the pores for photosynthetic gas ex- 
change (Fiscus et al. 2005; Roelfsema and Hedrich 2005). 
Adverse effects on plant photosynthesis were identified as 
a major factor limiting growth and yields of crop under 
high 0 3 concentrations. Another mechanism, which 
affects plants, is through the oxidative damage. The break- 
down of 0 3 in the apoplast is thought to lead to the 
formation of reactive oxygen species (ROS), such as super- 
oxide, hydrogen peroxide, and hydroxyl radicals. Such an 
0 3 -induced oxidative burst results in tissue damage that 
produces visible leaf damage leading to ROS-induced cell 
death process (Fiscus et al. 2005; Rao and Davis 2001; Frei 
et al. 2008; Pell et al. 1997; Black et al. 2007). 
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Studies found that the acute O3 concentration of 
560 nmol mol' 1 h could result in visible injury. This was 
obtained from accumulated exposure over duration of 7 
hours for eight consecutive days (Karlsson et al. 2004). 
Most of the researches carried out so far have been done 
with main focus on investigating the damages in bio- 
chemical, cellular, physiological, and morphological level 
of plants under chronic O3 exposure conditions for a 
relatively long period of a few weeks or few months 
(Cosgrove 2000; Zheng et al. 2000; Guidi et al. 2007; Gill 
and Tuteja 2010). 

As a way, to assess the effect of O3 on plant growth 
fast, we propose the use of Statistical Interferometric 
Technique (SIT), a non-invasive and non-contact optical 
interferometric method (Kadono et al. 2001; Kadono 
and Toyooka 1991). For measuring growth at higher res- 
olutions, attempts were made to use optical interferom- 
eter (Fox and Puffer 1976; Fox and Puffer 1977; Briers, 
1977; Jiang and Staude 1989). An optical interferometer 
involves the interference of reflected light from an optic- 
ally flat reference mirror and that from a plant. Usage of 
an interferometer was limited because of two main fac- 
tors. One is the complexity of the implementation and 
the other being, the optical property of the plant itself. 
Plant surface is highly scattering from both stationary 
and moving organelles on the surface and those within 
the leaf. The unwanted scattered light leads to the for- 
mation of a random pattern, called speckles (Dainty 
1984), obstructing the interferometer to achieve the ex- 
pected accuracy. In SIT, the totally random wave or 
speckle was used as reference instead of a plane wave 
generated from an optically flat mirror used in a conven- 
tional interferometer (Hariharan 1985). 

We demonstrated SIT to be capable of measuring 
plant growth at a spatial scale of nanometer (nm) and at 
a temporal scale of second (sec) under various environ- 
mental changes such as under O3 stress (Rathnayake 
et al. 2007) of Japanese red pine roots (Pinus densiflora 
Seibold & Zucc.) and infection of roots by two ectomy- 
corrhizal fungi (Rathnayake et al. 2008). Kobayashi and 
Kadono (2010) reported the effect of illumination condi- 
tions on the leaf-growth of rice crops by using SIT. 

In our recent investigation on the measurement of 
leaf-elongation of two different rice crops, Koshihikari 
and Fusaotome rice cultivars, it was revealed that the 
elongation rate contains random fluctuations at the 
nanometric spatial scale, nanometric intrinsic fluctua- 
tions (NIF). From characterizing these NIF by standard 
deviation, we found that NIF was influenced by envi- 
ronment such as the presence of O3. The NIF discovered 
differed between the cultivars and were also found to be 
present in different plants (Thilakarathne BLS, Rajagopalan 
UM, Kadono H, Yonekura T, High speed and high pre- 
cision optical interferometric technique to investigate 



instantaneous growth related changes of plant leaves, 
submitted). 

In this study, we focused on NIF and its dependence 
on repetitive cumulative O3 exposures of 7 hours for 
three consecutive days. O3 concentrations of 0 nl F 1 , 
120 nl 1, and 240 nl V were used. A Japanese rice culti- 
var, Koshihikari was chosen for its high sensitivity to O3 
(Yonekura 2011; Thilakarathne BLS, Rajagopalan UM, 
Kadono H, Yonekura T, High speed and high precision 
optical interferometric technique to investigate instantan- 
eous growth related changes of plant leaves, submitted). 

Results and discussion 

Leaf elongation rate under repetitive ozone exposures 

The elongation data of Koshihikari leaf obtained under dif- 
ferent time periods of minute to seconds under control O3 
concentration of 0 nl 1" are shown (Figure 1(A)-(C)). Over 
an hour, the elongation was in the order of tens of um and 
almost monotonically and slowly varying (Figure 1(A)). 
When confined to within less than a minute, random 
elongation fluctuations were observed in growth of a few 
hundred nanomters (Figure 1(B)). Under a few seconds as 
shown here for 6 sec, the elongation was almost linear and 
can be fitted to a line by linear regression (Figure 1(C)) to 
obtain the elongation rate. This elongation rate was nor- 
malized by the beam distance to calculate the relative 
elongation rate (RER). 

Figure 1(D) shows the relative elongation rate variation 
over 1 hour. The elongation rate was rapidly varying 
appearing as random noise fluctuations. It was found 
that the standard deviation (SD) of the fluctuations for 
the metal plate, which determine the noise level of the 
measurement system, was more than twenty times 
smaller than that for the leaf. 

Figure 2 shows the elongation rate fluctuations observed 
under three different O3 exposure conditions 0 nl 1" (top 
row), 120 nl I' (middle row), and 240 nl I' (bottom row) 
with each row containing data obtained on day 1 (left 
most), day 2 (middle), and day 3 (right most). The periods 
'Before) 'Within) and After' correspond to charcoal filtered 
(CF) air before O3 exposure, within the exposure, and CF 
air after the exposure, respectively, and are shown by sep- 
arating with broken lines. The fluctuation amplitudes de- 
crease with increase in concentration as well as due to 
cumulative exposures. 

In order to characterize these changes with exposure 
and repetition, the SD of these fluctuations was used as a 
norm and was referred as nanometric intrinsic fluctua- 
tions (NIF). Moreover, in order to reduce the common 
variations during a single session for a sample as well 
between different samples, normalization was done. The 
standard deviation of the fluctuations calculated under the 
initial one hour of CF air, SDj, e f ore , was used as the normal- 
izing factor. SD exposure and SD a fi e „ standard deviations of 
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Figure 1 Elongation of the leaf were continuously measured over 7 hours and shown at different time scales. Elongation of the leaf 
were continuously measured over one hour and shown at different time scales (A-C), and relative elongation rate for the leaf (D). The results 
were shown for an hour (A), a minute (B), and a few seconds (C). In (C), the straight line corresponded to a linear regression done to obtain the 
elongation rate. Elongation rate variation over an hour showed random fluctuations at nm spatial scales in conjunction with the monotonous 
macroscopic growth (D). Noise level of the system was measured with a metal plate as a sample, and the standard deviation of the fluctuations 
obtained with the metal plate was less than 0.1 nm mm" 1 sec" . 



the fluctuations under three hours O3 exposure, and under 
the last three hours of CF air exposure were normalized by 
SDbefore- The normalized fluctuations were called as nor- 
malized nanometric intrinsic fluctuations (NNIF). 

Normalized nanometric fluctuations for characterizing 
ozone effect 

NNIFs were calculated for all six samples under each O3 
concentrations of 0 nl 1"\ 120 nl 1"\ and 240 nl T 1 
and then averaged to give averaged NNIF. Figure 3 
showed the normalization procedure. The elongation 
rate fluctuations obtained from a leaf for three days 
are shown in Figure 3(A-C). The SDs during each 
time period of CF air before, 120 nl l" 1 O3 exposure, 
and CF air after were calculated and shown in Figure 3(D). 
Figure 3(E) showed the normalized standard deviations 
or NNIFs. 

Figure 4 showed the variation of NNIF as a function of 
repetition day under three different O3 conditions of 
0 nl 1" , 120 nl I" , and 240 nl I" . It can be clearly seen 
that, for 0 nl T 1 , there was almost no variation across 
the measurement duration and also across days. This 
indicated the measurement was relatively stable during 
the total period of three days. 



Results of repetitive ozone exposure for three days 

Both under 120 nl l" 1 and 240 nl f 1 , on the first day, a 
clear decrease in average NNIF could be seen within 
exposure period. On the first day of the exposure of 
120 nl 1" , measurement showed that there was no 
recovery to the pre-exposure level within the first day. 
However, on day 2, an increase in average NNIF could 
be seen at the start of the experiment, indicating a pos- 
sible recovery of the plant system from O3 stress on the 
first day. Again, with three hour exposure 120 nl f 1 O3 
on second day, there was a decrease in average NNIF 
but the decrease was getting smaller possibly due to 
acclimation of the plant to the O3 induced stress. On 
third day, there was actually increase again indicating 
the recovery with a reduced extent. 

For the case of 240 nl T 1 , the decrease in fluctuation 
amplitude was large, and with consecutive day exposures 
there was little difference in the decreases between the 
days. From Figure 4, we can see that there was a definite 
change in the amplitudes with repetitive exposure. We 
can also see some adaptation and the adaptation behav- 
ior depends on the concentration of O3. While for 
120 nl 1" , the plant's vitality is trying to get close to the 
pre-exposure level, the O3 effect becomes large enough 
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Figure 2 Results of nanometric intrinsic fluctuations (NIF). Results of nanometric intrinsic fluctuations (NIF) showing the effect of 0 3 under 
expeosure conditions of CF air (top row, (A-C)), 120 nl I" (middle row, (D-F)), and 240 nl F (bottom row, (G-l)) of 0 3 obtained on day 1 (leftmost 
column, (A, D, G)), day 2 (middle column, (B, E, H)) and day 3 (rightmost column, (C, F, I)). In figure, the durations 'Before', 'Within', and 'After' 
correspond to CF air before 03 exposure, within the exposure, and CF air after the exposure, respectively, and are shown by separating with 
broken lines. 
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Figure 3 Normalization procedure. Normalization process during each time period (before, within, and after) in three days for Koshihikari rice 
sample. Elongation rate variations shown as a function of time over three days (A-C). Variation of standard deviation as a function of time 
obtained during different time intervals (D). Variation of normalized standard deviations of elongation rate (E). Here standard deviation obtained 
at 'Before' on the first day was taken as the normalization factor. In all figures, (A)-(C), the durations 'Before', 'Within' and 'After' are shown by 
separating with broken lines. 
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Figure 4 Results of averaged standard deviation of NNIF. Results of averaged standard deviation of NNIF for Koshihikari rice cultivar showing 
the effect of 0 3 under three experimental conditions of 0 nl I" 1 , 120 nl I" 1 , and 240 nl f 1 0 3 exposures for three consecutive days. The error bars 
correspond to the SEs of the mean. Shown are means ± SE. The durations 'Before', Within' and 'After' are shown by separating with broken lines. 
For each treatment (CF, 120 nl F , and 240 nl F ) six rice samples from Koshihikari cultivar was used. All together 18 samples were investigated. 



that there was almost no change in the amplitudes with 
consecutive exposures for 240 nl 1" . 

Statistical analysis showed a significant difference 
(p < 0.01) between control and 120 nl 1" exposure, and 
control and after 120 nl 1" exposure. This was true for 
240 nl T too. Table 1 shows the decrease in amplitudes in 
percentages for three different conditions obtained under 
0 nl l' 1 , 120 nl F 1 , and 240 nl F 1 . 

NNIF and repetitive ozone exposure 

Due to the potential that SIT can monitor the changes 
not only before and after but also during the exposure of 
O3 making it a promising technique for monitoring 
continual changes in comparison to other techniques, 
we conducted an investigation with a continual expos- 
ure. In this study, leaf of Koshihikari cultivar was chosen 
for its relatively high sensitivity to repetitive exposure of 
O3 for three consecutive days. We found that there was 
no decrease in the normalized factor NNIF under con- 
trol, i.e., under an exposure of concentration 0 nl F 1 on 

Table 1 Reduction percentages in NNIF 



all three days. This suggests that our measurements were 
done in an almost robust condition during the total 
experimental period in relatively stable system and were 
not affected by external fluctuations like temperature 
and humidity. 

Under O3 exposure of concentration 120 nl F , there 
was respective immediate decrease in NNIF by 38% on 
the first day with no recovery back to pre-exposure state 
within the day. However, on the start of the second day, 
there was a recovery by 20% which started to show 
immediate decrease with exposure reaching a total de- 
crease to 67% on the second day. On the third day, for 
120 nl F recovery could be seen by 16% that remained 
almost constant throughout third day. The overall de- 
crease compared to pre-exposure was less than 25%. 
This indicated the recovery of the intrinsic adaptation of 
the system to continual exposure, although the recovery 
could never reach the pre-exposure state. 

In contrast, under 240 nl F , the decrease on the first 
day itself was large by 51% and there was no significant 



Repetition 
day 


% decrease in amplitude for 0 nl I" 1 


% decrease in amplitude for 120 nl I" 1 


% decrease in amplitude for 240 nl I" 1 


Before Within 


After 


Before 
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After 


Before 
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Day 1 


97 


96 




62 


64 




52 49 


Day 2 


99 101 


96 


84 


65 


67 


61 


50 45 


Day 3 


98 96 


95 


83 


75 


80 


4/ 


40 38 



Decrease in NNIF with repetitive three day exposures under 0 3 concentration of 0 nl l"\ 120 nl l"\ and 240 nl I" 1 . The decreases were those obtained during three 
time periods of before, within and after of 0 3 exposure. 
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change in decrease of NNIF between second and third 
days. While at the end of the second day the fluctuations 
decreased to 45% and on the end of the third day the 
fluctuations decreased to 38%. Therefore, the fluctua- 
tions became small that the deviations were much below 
the pre-exposure value at the first day and could never 
approach the pre-exposure level of NNIF. In actual case 
for 240 nl 1" , foliar damage could be seen only after a 
week of the three consecutive day exposures. The 
current investigation demonstrates the potential of SIT 
in immediate assessment of O3 exposure damage. 

Origin of NIF 

The observation using SIT is basically a macroscopic 
one since the probing region is the order of mm. The 
subnanometer accuracy of SIT makes it possible to 
monitor molecular or cell expansion changes as an 
integration of expansion changes over many cells. We 
hypothesize that the nanometric fluctuations are from 
possibly physiological changes accompanying plant 
growth. For growth of plant cells to happen, the cells must 
physically expand their restraining walls. At the same time 
the wall must preserve its mechanical integrity in the pres- 
ence of high turgor pressures. Spatial expansion may vary 
depending on the part of the plant and also on the various 
cell specific patterns of enlargement and divisions. Exist- 
ence of a variety of cell growth patterns and cell-specific 
enlargement leads to a variety of cell shapes. The spatial 
patterning is regulated by a variety of enzymatic processes 
that lead to modifications in the cell wall (Cosgrove 2000; 
Pantin et al. 2011). 

Sensitivity of NIF and 0 3 concentration 

Previous studies do indicate that the exposure to O3 
results in the reduction of both live regulating metabolic 
and hydraulic processes (Roelfsema and Hedrich 2005; 
Rao and Davis 2001; Black et al. 2007; Ruts et al. 2012; 
Wilkinson and Davies 2010; Kangasjarvi et al. 2005). As 
the nanometric fluctuations based on the results are 
hypothesized to be related to life regulating processes, 
we expected a decrease in NIF with increasing O3 accu- 
mulation. In SIT, the measured NIF decreased immedi- 
ately for large O3 concentration. 

In the current investigation, it was found that with 
repetitive exposure for three days of O3 under relatively 
small dosage of 120 nl 1" , it was possible to see the 
recovery while under a relatively large dosage of 
240 nl 1" for three days, the decrease in NIF was large 
and it was possible to observe some recovery on the 
second day that became almost null on third day. The 
irreversible damage in the leaf could be seen as foliar 
damage (Figure 5) only after a week. This agrees again 
with the earlier agreement of the results in reduction with 
gas exchange measures under 240 nl l" 1 (Thilakarathne 



BLS, Rajagopalan UM, Kadono H, Yonekura T, High speed 
and high precision optical interferometric technique to 
investigate instantaneous growth related changes of plant 
leaves, submitted). 

The physiological basis behind the origin of the fluctua- 
tions and the influence of the environment needs a further 
detailed investigation. The current study only focused on 
reporting the presence of such fluctuations and their 
sensitivities. 

Conclusions 

In this paper, we applied SIT for measuring the RER of a 
relatively low O3 tolerant Japanese cultivar, Koshihikari, 
to investigate its dynamic response to three hour O3 
exposure for three consecutive days at O3 concentrations 
of 0 nl I" 1 , 120 nl l' 1 , and 240 nl l' 1 . RER measurement at a 
temporal resolution of 0.5 sec and subnanometor accuracy 
revealed nanometric intrinsic fluctuations, NIF. To com- 
pare NIF across three days as well between different O3 
exposure concentrations, SD was used as a norm and 
normalized with those obtained on 1 st day before the 
exposure and averaged across six samples to give NNIF. 
NNIF was found to vary with the exposure concentrations 
of O3, and also NNIF to repetitive exposures revealed the 
dynamic response of possible recovery of the leaf under 
120 nl l" 1 O3 exposures. On the other hand, for an expos- 
ure of O3 concentration, 240 nl 1" was found to produce a 
large damage on the first day of the exposure. Although 
the following days' exposures indicated the evidence for 
the existence of recovery mechanisms, there was continual 
damage finally resulting in the appearance of foliar dam- 
age after a week. The possibility of monitoring ongoing 
response of the plant to environmental changes makes 
SIT an attractive tool. Conventional techniques such as 
gas exchange measures basically obtain data before and 
after change in environment and it is not possible to 
measure in a dynamic way. 

Methods 

Experimental system of SIT 

Figure 6(A) shows a schematic of the statistical interfer- 
ometer assembled on an optical bench (Meiritsu seiki, 
Japan) used to measure plant growth. A photograph of 
the optical system is given in Figure 6(B). The light 
emerging from a He-Ne laser of wavelength of 633 nm 
(GLG 5400, NEC corporation, Japan) was first passed 
through a neutral density filter (NDF) (F71N-2, Suruga 
Seiki, Japan) to adjust the light intensity. The light was 
further divided into two beams by a specially designed 
prism PI (Suruga Seiki, Japan). Two beams normally 
illuminated two points on the surface of a growing leaf. 
In the experiment, third or fourth leaf of the rice plant 
mounted in the chamber with a custom-made holder 
was used (Figure 6(C)). Precaution was taken to avoid 
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Figure 5 Foliar damage observed after a week. Foliar damage observed after a week for three days with 0 3 exposure of 0 nl I" 1 (A) and 
240 nl I" (B). The yellowish spot seen on illuminated spots of the leaf and on the tip resulted from the consecutive exposure of 0 3 for three 
consecutive days. Black color circle showed foliar damage area. The damage appeared after a week from the start of the exposure. The damage 
was visible only under 240 nl I" 1 and not under 120 nl I" 1 although under both levels the damage was irreversible. The illumination spots were 
identified based on the flour remaining on the leaf after finishing the exposure experiments. 



damage to the leaf by having soft material such as cotton 
to cover the surface of the holder. 

On illuminating two points on the leaf, two independent 
speckle fields were generated that interfere to make a ran- 
dom interference speckle pattern as shown in Figure 6(A). 
The speckle patterns were acquired by a CCD camera 
(XC-75, Sony Corp., Japan) set at a distance of 83 mm 
from the leaf. An interference filter (IF) was attached in 
front of the CCD camera to suppress the background light 
and select the probing light. The IF worked as a band-pass 
filter having a central wavelength of 633 nm with a band- 
width of 10 nm. A polarization filter (PF) was also used to 



maximize the interference signal. The distance d between 
the points of illumination on the leaf was set to be 3 mm. 
d can be adjusted by moving the prism P2 (Suruga Seiki, 
Japan) along the optical axis. P2 was mounted on a piezo- 
electric transducer (PZT) stage (E-620, Physik Instrumente, 
Germany) was used to give a certain phase shift between 
the interfering light beams. A frame grabber (VEC-PRO, 
IMPERX, USA) was used to acquire and store the interfer- 
ence speckle patterns every 500 msec. 

Further, the probing area was covered with wheat flour 
to avoid light penetrating inside the leaf. Such penetration 
of light would result in generation of dynamic speckles 
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Figure 6 A schematic of the experimental system of SIT. A schematic of the experimental system of SIT (A) used for measuring plant growth 
placed in an 0 3 chamber with a photo of the real system on an optical bench (B). Third or fourth leaf was stably fixed using a sample holder (C) 
with precaution taken to avoid any damage to the leaf. The leaf length was about 1 5 to 20 cm. The probing points on the leaf were set to be 
around 4 cm from the apex of the leaf. The distance d between the points of illumination on the leaf was set to be 3 mm, and d can be 
adjusted by moving the prism P2 along the optical axis. The CCD camera was placed at a distance of 83 mm from the probing area, and the 
speckle interference patterns (A) were continuously recorded. Here, the notations indicate, NDF: neutral density filter, PZT: piezoelectric 
transducer, PI and P2: prisms to generate two parallel probing beams, 0: angle between the probing and observing directions. 



from moving organelles within the leaf degrading the 
measurement accuracy of SIT. Such a dynamic speckle 
pattern is known as biospeckle (Aizu and Asakura 1996). 
Three halogen lamps with fiber optic light guides 
(PHL-150, Mejiro Precision Ltd., Japan) were used to pro- 
vide the light intensity of 500 umolm" sec" . 

SIT developed in our lab. is now available as a commer- 
cial product from Toyo Seiki Seisaku-sho, Ltd., Japan. 

Principle of SIT 

SIT is a non-contact optical interferometric technique 
for measuring displacements at an accuracy of subnan- 
ometer and a temporal resolution of second (Kadono 



et al. 2001; Kadono and Toyooka 1991). Principle of stat- 
istical interferometry is different from that of conven- 
tional interferometric techniques. In a conventional 
interferometer, the light reflected from an optically flat 
reference mirror and reflected light from the sample 
under observation interfere to give an interference pat- 
tern. The phase change is derived from the interference 
pattern. It corresponds to the change of the optical path 
such as a displacement of the sample. In order to 
achieve subnanometer accuracy, stringent conditions 
were imposed on the interferometer that all of the 
uncertain factors in the interferometer, e.g. distortion of 
the reference mirror, vibration of the optical components, 
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air turbulence, electric noise of the circuit and so on, 
has to be eliminated. Especially on the reference arm, 
it is required to make the resulting reference wave- 
front with no uncertainty such as deformity of refer- 
ence mirror. 

In contrast, in statistical interferometer, a completely 
random wavefront called a fully developed speckle field 
is used as reference and is generated when an optically 
rough surface is illuminated by a laser light. Phase of the 
speckle field has a uniform probability density function 
(PDF) that takes a constant probability density of 1/2it 
over the range from - re to n Characteristics of such 
speckle phase is remarkably stable and reliable and can 
play the role of a standard. To understand the principle 
of the method, one has to know the fundamental prop- 
erties of the speckle field and its statistics that are given 
in (Additional file 1). 

For every 500 milliseconds, the speckle interference 
patterns were acquired by a CCD camera and stored in a 
frame memory. Basically, when there is an in-plane 
elongation of Ax between the two illuminating points of 
the object, the optical path difference AL between the 
two interfering speckle fields can be expressed by, 

AL = Axs'm8, (1) 

where 6 is the angle between the illumination and the 
observation directions. Due to the elongation Ax, the 
phase difference between the two interfering light fields 
changes by W and is given in terms of AL as: 

V = (2n/X)AL. (2) 

SIT data analysis to obtain elongation rate 

In our SIT algorithm, for deriving the unknown object 
phase the statistical property of the speckle phase, i.e., 
the total randomness of the speckle phase was used as a 
constraint. A detailed basic SIT algorithm, and its exten- 
sion developed by the authors, and derivation of the 
unknown object phase are given in (Additional file 1). 
The whole processing could be performed to get the 
true object phase with very high accuracy without any 
time consuming iterative procedure. 

From the calculated phase change, corresponding 
elongation was calculated in real time. From the derived 
elongation over a period r of 5.5 sec, the slope of the 
elongation data was obtained by linear regression. Rela- 
tive elongation rate (RER) (nm mm' sec' ) is defined by, 

Relative Elongation Rate = Elongation / rd . (3) 

The following procedure was applied to obtain the 
RER and its fluctuations under control and O3 expo- 
sures: (a) Linear regression was done over 5.5 sec to obtain 
slope or to obtain instantaneous RER; (b) Repetition of 



procedure (a) for every 5.5 sec across the whole session of 
7 hours; (c) The calculated RER was passed through a 
moving average filter with window size 165 sec to ob- 
tain smoothened data; (d) The smoothened data was 
subtracted out from the RER data of step (b) to elim- 
inate a long-term trend. All calculations were done 
with Matlab (R2011b, MathWorks, USA). Statistical 
analyses were performed with the SPSS (Ver.19) stat- 
istical package. To analyze the statistical significance 
of the standard deviations obtained under different 
concentrations of O3, i-test was used. 

Plant cultivation and 0 3 exposure 

This research was carried out from December 2010 to 
May 2011 in Saitama University, Japan. Thirty to forty 
days old Koshihikari rice cultivar was used. We chose 
the age based on earlier studies (Inada et al. 2008; Yama- 
guchi et al. 2008). That study showed that during the 
vegetative period of thirty to forty days, exposing the 
cultivar to O3 had significant reduction on the relative 
growth rate, net assimilation rate, and photosynthetic 
rate. The culturing conditions were as follow: 

1. The seeds of Koshihikari rice plant were 
pre-germinated in dark environment 
for four days. 

2. Seedlings were transferred to separate plastic 
containers (600 ml volume cups) with two plants 
per container. Plastic containers were filled with the 
same soil sample with the required fertilizer, up to 
300 ml of the cup. 

3. The seedlings were placed in a growth chamber 
(Conviron, Controlled Environmental Ltd, 
Winnipeg, Manitoba, Canada.) for around three 
weeks. The conditions of the chamber were 
made to follow a day/night cycle of: 

a) 12 hours/12 hours; 

b) Light intensities of 260-350 umol m' 2 sec' 1 / 
0 umol m' 2 sec' 1 ; 

c) Temperatures were 25°C/18°C; 

d) Relative humidity was kept at 55%-65% 
throughout a day; 

Figure 7 shows a schematic of O3 exposure system 
consisting of O3 chamber (33.5 x 33.5 x 58 cm 3 ) made 
up of glass and metal, an air pump, charcoal filters, an 
O3 generator (OES 10A, Dylec Inc.), and an O3 monitor 
(Model 1150, Dylec Inc.). The air from the air pump was 
passed through a charcoal filter (CF) to remove the 
ambient O3. The same inlet and outlet in the cham- 
ber were used for O3 and the charcoal filtered air. O3 
concentrations were monitored with the O3 monitor, 
and the signal was fed back to the O3 generator to 
control O3 concentration precisely. The O3 feedback 
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Figure 7 A schematic of experimental system. A schematic of experimental system for introducing Charcoal filtered air (CF air) and 0 3 . Signal 
from the 0 3 monitor was fed back to the 0 3 generator to control precisely the concentration of 0 3 inside the chamber. Charcoal filter CI was 
used to remove the ambient 0 3 in the environment. C2 and C3 were used to remove all 0 3 before releasing the air into the room. Flow rate was 
10 l/min. Three halogen lamps with fiber optic light guides were used to provide light intensities of 500 umolirf 2 s" 1 from the top of 0 3 chamber 
(not shown explicitly). 



control system had an accuracy of around 0.5% during the 
experiment. 

In O3 exposure measurements, the O3 generation 
system approximately took 10 minutes and 14 mi- 
nutes to reach the set O3 concentrations, 120 nl l" 1 
and 240 nl 1" , respectively. Similarly the control sys- 
tem took approximately 10 minutes and 14 minutes to 
bring O3 concentration to 0 nl I" 1 from 120 nl I" and 
240 nl 1" , respectively. O3 concentrations used here 
were chosen based on the common criteria used in 
Japan. The condition of 120 nl 1" corresponds to issuing 
warning while 240 nl T 1 corresponds to issuing a serious 
warning. 

0 3 exposure timing protocol 

At first, continuous CF air (10 liters/min.) was fed 
into the chamber for two hours. During this initial 2 
hours, the stabilization of the whole system was care- 
fully monitored, and then the data acquisition started. 
All the experiments were carried out at fixed times 
from 11.00 a.m. to 6.00 p.m. The entire experimental 
duration took 7 hours. The timing protocol used was 
as follows: 

1. CF air was provided for the first 1 
hour (control); 

2. Next 3 hours O3 was introduced into 
the chamber; 

3. For the last 3 hours, again CF air was used 
to fill the chamber. 



Three conditions namely, the control with O3 concen- 
tration being 0 nl T 1 , O3 exposure conditions of 120 nl 1" 
and 240 nl l' 1 , were used. For each O3 exposure condition 
of control, 120 nl T 1 , and 240 nl T 1 , six samples were used. 
In control experiments, only CF air was provided during 
the entire duration of 7 hours. 

Both control and O3 exposure experiments were con- 
ducted for three days consecutively. At the end of the 
each day, the plant used for experiment was kept in the 
growth chamber to be used on the following day. On 
second and third days, the same probing area of the leaf 
as of the first day (identified through the wheat flour on 
the leaf surface), was used to get the SIT measurements. 
The accumulated exposures over the threshold 40 nl 1" 
(AOT 40) for three consecutive days were 240 nl 1" .h, 
480 nl r\h, and 720 nl l'\h for O3 exposure of 120 nl I 1 , 
and 600 nl l^.h, 1200 nl l^.h, and 1800 nl I~\h for O3 
exposure of 240 nl I" 1 . 

Additional file 



Additional file 1: SI. Fundamental statistical properties of the speckle 
field and the effect of non-spatial uniformity of object. S2. Algorithm of 
Statistical Interferometry. S3. Expansion of dynamic range in statistical 
interferometry. 
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SIT: Statistical interferometric technique; PDF: Probability density function; 
CCD: Charge coupled device; CF: Charcoal filtered; NIF: Nanometric intrinsic 
fluctuations; NNIF: Normalized nanometric intrinsic fluctuations; SD: Standard 
deviation; PZT: Piezoelectric transducer. 
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